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Повышение качества 

и чувствительности анализа 
низкоконтрастных изображений методом 
виртуальной оптической эллипсометрии 


Рассмотрены информационные возможности нового метода повышения качества и чувствительности 
сегментации низкоконтрастных изображений. Суть метода основана на переходе от амплитудной 
(яркостной) модуляции к фазовой (угловой) модуляции в комплексной области, что открывает 
возможность обобщения идеи оптической эллипсометрии на область цифровых методов обработки 
изображений. Представлены результаты численного моделирования и реальной проверки 
работоспособности метода. 


Введение 


При разработке новых методов компьютерного зрения допустимо использование 
виртуальных аналогов физических методов обработки радиофизических и оптических 
сигналов и полей и их математических моделей. Известно, что наибольшую чувстви- 
тельность в радиофизических и радиотехнических измерениях обеспечивают резонансные 
методы, а в оптических — голографические, интерференционные и эллипсометрические. 

Все указанные измерительные методы базируются на использовании волновых 
и поляризационных свойств когерентного электромагнитного излучения. Для возмож- 
ности использования данного подхода необходимо преобразование исходного изображе- 
ния к такому виду, который бы обеспечивал возможность использование идей и 
математического формализма указанных физических методов. 

В работе [1], [2] нами было предложено использование модуляционного преоб- 
разования вида 


‚ (ху) 


Ву) =КЦьуе” * =Ва(ь у)}+ Лифксь у} = Во у)е®”, (1) 


где Г(х, у) — исходное изображение, а А — модуляционный параметр. 


Использование модуляционного преобразования (1), обеспечивающего переход 
в пространство комплексных яркостей, открывает еще одну дополнительную возможность 
анализа изображений, основанную на проведении аналогий с поляризационными свой- 
ствами волновых полей. Как это хорошо известно, наиболее чувствительным оптическим 
методом исследования пленочных структур является метод оптической эллипсометрии, 
теоретические и практические аспекты которого подробно изложены в [3], [4]. 


«Штучний 1нтелект» 32008 233 


Ахметшина Л.Г. 


Цель данной статьи — демонстрация возможности повышения качества и раз- 
решающей способности визуального анализа низкоконтрастных цифровых изображений 
на основе проведения формальной аналогии с физическим методом оптической эллипсо- 
метрии. 

Постановка задачи. Исследование информационных возможностей метода 
виртуальной оптической эллипсометрии применительно к трем классам изображений: 
обычное яркостное изображение, цветное изображение (сюда относится и случай 
цветового КОВ кодирование изображений); мультиспектральное изображение (число 
диапазонов/параметров больше трех). 


Виртуальная эллипсометрия однопараметровых 
изображений 


Рассматривая результат модуляционного преобразования А(х, у) в виде аналога 


двумерной проекции трехмерного векторного поля В, можно естественным способом 
ввести и векторное поле его ротора [5], [6] 


то В= 


ОК Е 
дк. ОК, н [№ @ ]= ‚› _ В, ©) 


9 08 2 0х 0х ду 


Выражение (2) является справедливым, если частные производные А(х, у) являются 


непрерывными, что, безусловно, выполняется для классов изображений, ЯВЛЯЮЩИХСЯ 
предметом исследования. Поскольку компонента и является неизвестной, ТО В ВЫ- 


ражении (2) может быть определена лишь третья компонента 7ой В. 


ОК, дЮ 
А. (х, у) = —. о (3) 


которую можно рассматривать в качестве дополнительного информативного параметра, 
ортогонального составляющим А, и К. 


Использование такого подхода позволяет сопоставить каждому пикселю исходного 
изображения три ортогональные компоненты, что открывает возможность использования 
волновых оптических преобразований, вытекающих как из самой сути модуляционного 
преобразования (А можно рассматривать в качестве виртуального аналога длины волны), 
так и из трехмерной ортогонализации пространства синтезированных информатив- 
ных признаков. 

Действительно, если рассматривать компоненты А, и К, в качестве соответствую- 


щих амплитудных составляющих эллиптически поляризованной ВОЛНЫ, ТО каждой 
точке плоскости «х-у» можно сопоставить поляризационный эллипс, главная ось 


которого образует азимутальный угол и с осью х: 


1 25,6, 
и= - ОЕ Е с050), (4) 


у 
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где Би Б, рассматриваются как модули компонент К, и К, соответственно. Основ- 


ная тонкость предлагаемого подхода заключается в способе вычисления фазового 
угла д=д,(х,у) между компонентами А, и К, , рассчитываемого как 


А. (2 
А 05) 


= (у =е” * , (5) 


где в качестве А. рассматривается нормированное значение компоненты ротора 


векторного поля В. Помимо этого, отношение меньшей полуоси поляризационного 
эллипса с к большей полуоси 4 можно выразить как функцию вспомогательного 
угла (угла эллиптичности) у, так что с/4 = 1х ‚ где 


1 ОЙ 
д= а то). (6) 


Таким образом, первоначальному изображению можно сопоставить два новых: 
азимутальное изображение Ч(х,у) и эллиптическое изображение у(х,у), которые 
характеризуются нелинейной комбинацией всех трех синтезированных ортогональ- 
ных параметров. 

Столь же полно, в рамках рассматриваемого подхода, анализируемое изображение 
может быть охарактеризовано параметрами Стокса (рис. 1), определяемыми как [7] 


55 ЕБЕ+Ь; 5, = 2,5, с036 = 5 032 узш2и ; 


В, 5; 53 =26.6, зтд = 5, зт2х. (7) 


Рисунок 1 — Сфера Пуанкаре 
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Параметр 5, пропорционален интенсивности, а 5,,5,,5, интерпретируются, как 
декартовы координаты точки на сфере радиусом 5,, известной как сфера Пуанкаре. 
Долгота и широта этой точки равны соответственно 2 и 2х. 

Исследования были проведены на многих образцах низкоконтрастных изобра- 
жений. Ниже представлены некоторые из полученных результатов в сравнении с 
известными преобразованиями (значение д =1). 

Маммографическое изображение. На рис. 2 представлено медицинское радио- 
логическое изображение и результаты его обработки известными методами [8], [9]. 


— = 
9 Гы 


а б В 
Рисунок 2 — Маммографическое изображение: а — оригинал; б — результат метода 
эквализации гистограмм; в — модуль градиентного отображения 


Из рассмотрения рис. 2 видно, что применение известных методов анализа низко- 
контрастных изображений не позволяет детализировать потенциальную область 
возможного интереса. Использование эллипсометрических характеристик (ху) и 


Я(х, у) (рис. 3) позволяет повысить детализацию визуального анализа. 


Ч (х, у) Хх, у) 


Рисунок 3 — Эллипсометрические характеристики 


Из рассмотрения рис. 3 следует, что синтез виртуальных эллипсометрических 
характеристик на основе комбинации методов модуляционного преобразования И 
теории виртуальных векторных полей позволяет повысить чувствительность И 
резкость однопараметровых (яркостных) низкоконтрастных изображений, при этом 
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анализ азимутальных синтезированных характеристик Ч(х,у) больше подходит для 
решения задач сегментации, тогда как синтез эллиптических характеристик х(х, у) 
позволяет повысить чувствительность и резкость исходных изображений. 

Инфракрасное изображение. На рис. 4а представлено инфракрасное изобра- 
жение работающей интегральной микросхемы, а на рис. 46 — результат применения 
метода эквализации гистограмм. 


и 


а б }: 


Рисунок 4 — Инфракрасное изображение: а — оригинал; б — результат метода 
эквализации гистограмм; в — эллиптическая характеристика у(х,у) 


Из рассмотрения рис. 4 следует, что метод эквализации гистограмм в данном 
случае является неэффективным, тогда как синтез эллиптической характеристики 
Х(х, у) позволил выявить все топологические особенности анализируемого изображения 


и, в частности, выявить характерную ступенчатость его структуры связанной с особен- 
ностями его формирования на основе использования ПЗС матрицы. 

Виртуальная эллипсометрия цветных изображений. В качестве цветного 
может рассматриваться любое изображение, представленное тремя параметрами 
= (%, у), 1, =1.(х,у) и 1, =1[\(х,у) (включая формат КСВ -— или Н$У кодирова- 
ние). Данное обстоятельство открывает возможность непосредственного использования 
метода виртуальной эллипсометрии для повышения чувствительности и разрешаю- 
щей способности анализа цветных изображений. 

Синтез трех изображений виртуальных углов эллиптичности 7, = (х,у), 
ДЕЛ (ху) и д, = Х,(х, У) возможен на последовательном использовании следующих 


комбинаций: 
д = ЖИ, 1,6. = ЛИ, 
д = д» Ь1ь,6, = ЛА,}; (8) 
дз = 1,1, = И, }, 


где д, вычисляются на основе использования выражения (5). 


Следует иметь в виду, что в отличие от предыдущего случая (однопараметро- 
вого изображения), в котором все три компоненты синтезированного векторного поля 
являлись ортогональными, в выражении (8) любая комбинация исходных изображений 


у м. ‚} этим свойством не обладает. Это важное обстоятельство не позволяет непо- 
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средственно использовать эллипсометрический метод. В выражении (8) необходимо 
использовать ортогонализированные аналоги | Г.„ | тогда как в выражениях 0, 
и 


от; > ом 
используются прямые зависимости 1 1 задающие величину фазового сдвига между 
ортогональными составляющими | ий, Ё для фиксированного значения параметра 
1 у 


Я в выражении (5). Ортогонализация комбинаций И Ш. ы производилась методом 
Грама — Шмидта [10]. 

Проверка работоспособности метода для реальных изображений в целях сокращения 
объема изложения представлена на примере использования только эллиптической 
характеристики 7(х, у) и одной комбинации { Пана о =) р поскольку резу- 


льтат является типичным и для других вариантов. 

На рис. 5 приводится цветной фрагмент изображения акватории Черного моря 
(рис. 5 а) и результат эллипсометрического синтеза, в сопоставлении с результатами 
преобразования известными методами. Чрезвычайно низкая контрастность исходного 
цветного изображения обусловлена тем обстоятельством, что перепады яркостей во 
всех трех диапазонах находятся в пределах 4-- 6 %. В данном случае методы нечетких 
С-средних и самоорганизующихся карт Кохонена обеспечивают болыпую степень 
детализации, по сравнению с исходным снимком, но они лишь сегментируют исходное 
изображение, однако не могут повысить степень пространственной разрешающей способ- 
ности. Если анализировать характеристику х(х, у) ‚то в окрестности, обозначенной «стрел- 


кой», явно выделяются детали, невидимые ни на одном из других изображений рис. 5. 


8. 


42 


а б р: Г 


Рисунок 5 — Сопоставление различных методов преобразования: а — оригинал 
(в градациях серого); 6- 7(х, у); в- нечеткие С-средние (10 кластеров); г- карта 


Кохонена 


Чтобы показать, что характеристика у(х,у) действительно обеспечивает повы- 


шение пространственной разрешающей способности, целесообразно рассмотреть 
пространственно-амплитудные срезы исходных трех диапазонов и синтезированной 
характеристики 7(х,у), представленных на рис. 6 (в данном случае представлены 
100-ые столбцы). Тот факт, что амплитудно-частотные характеристики формирую- 
щих фильтров усиливают высокие пространственные частоты, свидетельствует об 
увеличении пространственной разрешающей способности при использовании эллип- 
сометрических характеристик. 
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Рисунок 6 — Пространственно-амплитудные срезы: а, 6, в — три диапазона исходного 
изображения; г- у(х, у) ; амплитудно-частотные характеристики формирующих 
фильтров: д- вход «а», выход «г»; е — вход «б», выход «г» 

На рис. 7 представлены результаты синтеза эллипсометрических характеристик 
в ситуации, полностью аналогичной вышеизложенной, но без использования предва- 
рительной процедуры ортогонализации составляющих /\(х,у)и [.(х, у). 


а | 


ЧФ (у) Я, У) 
Рисунок 7 — Эллипсометрические характеристики без использования 
ортогонализации составляющих /| и /), 

Результаты, представленные на рис. 7, свидетельствуют, что без использования 
предварительной процедуры ортогонализации составляющих (в данном случае / и /,), 
повышения пространственной разрешающей способности (и как следствие, чувстви- 
тельности визуального анализа) не происходит. Поэтому процедура ортогонализации 
является необходимой в рамках данного метода. 

Эллиптическая характеристика 7(х, у) обладает еще одной интересной особен- 
ностью, сближающей ее в некоторой степени с методом эквализации гистограмм. На рис. 8 
представлено аэрокосмическое КСВ - изображение и результат применения к каждому 
диапазону метода эквализации гистограммы и аналогичные его эллиптическая характе- 


ристика у(х,у). 
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Применение метода эквализации гистограмм улучшает яркость и контрастность 
цветовой палитры, по сравнению с исходным изображением, что приводит к повышению 
чувствительности при визуальном анализе. Применение нового метода (эллиптической 
характеристики 7(х, у) ) позволяет, кроме перечисленных факторов, увеличить также 
пространственную разрешающую способность, что выражается в детализации изоб- 
ражения. 


Рисунок 8 — ВОВ изображения, конвертированные к уровням серого: 
а — оригинал; 6 — метод эквализации гистограммы; в — характеристика у(х, у) 


На рис. 9 представлены гистограммы распределения яркостей трех диапазонов 
исходного изображения (рис. 8а), а на рис. 10 гистограммы распределения яркостей 
трех синтезированных эллиптических изображений (рис. 8в). 
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Рисунок 9 — Гистограммы распределения яркостей трех диапазонов исходного 
изображения (рис. 8а) 
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Рисунок 10 — Гистограммы распределения яркостей трех диапазонов для 
синтезированного эллиптического изображения (рис. 8в) 


Видно, что последние распределены по всему диапазону градаций яркостей и 
весьма близки, за исключением 7/,(х, у), к равномерному распределению, что и обус- 
ловливает лучшую степень зрительного восприятия результатов КОВ кодирования. 
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Виртуальная эллипсометрия мультиспектральных изображений. Как пример 
подобного типа изображений на рис. 11 показаны первые шесть изображений 8-мер- 
ного мультиспектрального ансамбля. 


. — 


Рисунок 11 — Первые шесть изображений мультиспектрального ансамбля 


Большое количество исходных диапазонов затрудняет синтез эллипсометрических 
характеристик для всего ансамбля и в этой связи возникает задача сжатия исходных 
данных, в наилучшем варианте до трех изображений. 

Для этой цели нами был использован метод сингулярного разложения, причем, 
степень сжатия исходного ансамбля определялась на основе анализа спектра нормиро- 
ванных сингулярных чисел, который показал, что первые три «собственных» изображения 
ансамбля содержат более 99 % информации мультиспектрального ансамбля (каждое 
«собственное» изображение представляет собой линейную комбинацию всех исходных 
изображений ансамбля, но с разными значениями коэффициентов), что дает основание 
для вывода о достаточности их использования для непосрественного цветового ВОВ-ко- 
дирования (рис. 12г). Из рассмотрения рис. 12 следует, что цветовое кодирование трех 
«собственных» изображений ухудшает визуальное восприятие полученного изображе- 
ния, по сравнению с «собственными» изображениями. 


а б в Г 


Рисунок 12 — Преобразование мультиспектрального ансамбля: а, 6, в — первые три 
«собственных» изображения; г — результат К ОВ-кодирования (в градациях серого) 
на их основе 


На рис. 13 представлены результаты синтеза трех эллиптических изображений и 
их цветового кодирования, полученные на основе использования трех «собственных» 
изображений, показанных на рис. 12а -— в. 

Сопоставление этих результатов (рис. 12г с рис. 13) показывает, что использование 
эллиптических характеристик первых трех «собственных» изображений позволило не 


«Штучний 1нтелект» 32008 241 


Ахметшина Л.Г. 


только применить эллипсометрическую методику к мультиспектральному ансамблю, но 
и обеспечить высокую степень детализации при КСВ-кодировании на основе трех 
синтезированных эллиптических характеристик. 


ре КСВ 


Рисунок 13 — Эллиптические характеристики и их цветовое ВКС В-кодирование 
(в градациях серого), полученные на основе использования первых трех 
«собственных» изображений исходного мультиспектрального ансамбля 


Заключение 


Метод виртуальной оптической эллипсометрии, в основе которого лежит модуля- 
ционное преобразование яркости изображений, открывает возможность проведения 
аналогий с поляризационными свойствами волновых полей с использованием соответст- 
вующего математического аппарата. Данный подход обеспечивает повышение чувст- 
вительности и разрешающей способности на основе нелинейного синтеза композиционных 
эллипсометрических характеристик: азимутального и эллиптического углов. 

Практическими особенностями его применения являются следующие факторы: 
— в случае обычного яркостного изображения, результат модуляционного преобразо- 
вания рассматривается в виде аналога двумерного векторного поля, а третья ортого- 
нальная компонента — в виде 2 -компоненты ротора этого векторного поля, что позволяет 
непосредственно синтезировать эллипсометрические характеристики; 

— при анализе цветных изображений (три измерительных канала), первым шагом 
является проведение предварительной процедуры ортогонализации двух из исходных 
изображений с применением модуляционного преобразования к третьему (процедура 
осуществляется в трех последовательных комбинациях), рассматриваемых в виде аналога 
двух поляризационных амплитуд и фазового угла между ними соответственно с после- 
дующим синтезом на их основе трех эллиптических характеристик и осуществления их 
цветового ВОВ-кодирования, при этом чувствительность, разрешающая способность 
и степень цветового зрительного восприятия превосходят характеристики оригинала; 
— для мультиспектрального ансамбля, когда число измерительных каналов больше трех, 
необходима процедура его сжатия до трех изображений, в целях использования методи- 
ки, применяемой для анализа цветных изображений; эффективным для сжатия является 
метод сингулярного разложения, причем в качестве трех основных («собственных») 
изображений принимаются те, которые соответствуют первым трем максимальным син- 


242 «Искусственный интеллект» 32008 


Повышение качества и чувствительности анализа низкоконтрастных изображений... 4 А 


гулярным числам — синтез эллиптических характеристик на основе этих «собственных» 
изображений, обеспечивает более высокую степень детализации результирующего цве- 
тового КОВ-кодирования по сравнению с прямым КСВ-кодированием трех «собст- 
венных» изображений. 
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ЛА. Ахметшина 

двищення якост! та чутливост! аналзу низькоконтрастних зображень методом в!ртуально! 
оптично! елепсометри 

Розглянуто 1нформащйн! можливост! нового методу шдвищення якост! 1 чутливост! сегментаци 
низькоконтрастних зображень. Суть методу полягае на переход! в1д амплтудно! (яскравоТ) модуляци 
до фазово! (кутовоТ) модуляцИ в комплексно! област!, що в1дкривае можлив!сть узагальнення 1дет 
оптично! елепсометрй на галузь цифрових методв обробки зображень. Представлен! результати 
реально! переврки працездатност! методу. 


Г.С. АКйтезйта 

Тпсгеазшо ОпаНбу апд Зепзвуйу Апа!уз5 о? Гоуу Сопбга$Е Гтасез Вазе4 оп Уйиа! Орйса! Ерзотейгу 
Метод 

шЕРгпайоп роз$ШЯез оЁ пе тео оЁ шсгеазше диаШу апа зеотеганоп зепзуйу оЁ 1о\\ сопгаз 
парез аге соп$14еге4. Тве таш 14еа оЁ Фе тео 1$ Базе оп Ше тапз оп Нот тагпаае (Биетез$) 
тодайоп, 0 рВазе (апоШаг) тодиаНоп ш а сотр/ех 5расе. [ ргоу!4ез а роз5Ииу г вепега|ханоп 
1еа оЁ орйса! еШрзотейу 10 дотат о# 42а] те о4$ оЁ итазе ргосеззте. Кези$ оЁ геа] {е5Нпо оЁ Те 
тео роз Шиез аге ргезещеч. 


Статья поступила в редакцию 05.06.2008. 
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